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Abstract 

We introduce and solve a model which considers two coupled networks 

growing simultaneously. The dynamics of the networks is governed by the 

new arrival of network elements (nodes) making preferential attachments to 

pre-existing nodes in both networks. The model segregates the links in the 

networks as intra-links, cross-links and mix-links. The corresponding degree 

distributions of these links are found to be power-laws with exponents having 

coupled parameters for intra- and cross- links. In the weak coupling case the 

model reduces to a simple citation network. As for the strong coupling, it 

mimics the mechanism of the web of human sexual contacts. 
PACS numbers: 02.50.cw, 05.40.-a, 89.75Hc. 
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I. INTRODUCTION 



Today, with a vast amount of publications being produced in every discipline of scientific 
research, it can be rather overwhelming to select a good quality work; that is enriched with 
original ideas and relevant to scientific community. More often this type of publications 
are discovered through the citation mechanism. It is believed that an estimate measure for 
scientific credibility of a paper is the number of citations that it receives, though this should 
not be taken too literally since some publications may have gone unnoticed or have been 
forgotten about over time. 

Knowledge of how many times their publications are cited can be seen as good feedback 
for the authors, which brings about an unspoken demand for the statistical analysis of 
citation data. One of the impressive empirical studies on citation distribution of scientific 
publications |E] showed that the distribution is a power-law form with exponent 7 ~ 3. The 
power-law behaviour in this complex system is a consequence of highly cited papers being 
more likely to acquire further citations. This was identified as a preferential attachment 
process in [0. 

The citation distribution of scientific publications is well studied and there exist a number 
of network models [Q-|5j to mimic its complex structure and empirical results to confirm 
predictions. However, they seem to concentrate on the total number of citations without 
giving information about the issuing publications. The scientific publications belonging to 
a particular research area do not restrict their references to that discipline only, they form 
bridges by comparing or confirming findings in other research fields. For instance most Small 
World Network Models []7|~fl presented in statistical mechanics, reference a sociometry article 
|T0| which presents the studies of Milgram on the small world problem. This is the type of 
process which we will investigate with a simple model that only considers two research areas 
and referencing within and across each other. The consideration of cross linking also makes 
the model applicable to the web of human sexual contacts fTl| : l^|, where the interactions 
between males and females can be thought of as two coupled growing networks. 
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This paper is organized as follows: In the proceeding section the model is defined and 
analyzed with a rate equation approach U^,|TJ]. In the final section discussions and com- 
parisons of findings with the existing data are presented. 

II. THE MODEL 

One can visualize the proposed model with the aid of Fig. ([!]) that attempts to illustrate 
the growth mechanism. We build the model by the following considerations. 

Initially, both networks A and B contains mo nodes with no cross-links between the nodes 
in the networks. At each time step two new nodes with no incoming links, one belonging to 
network A and the other to B, are introduced simultaneously. The new node joining to A 
with rriA < m outgoing links, attaches paa fraction of its links to pre-existing nodes in A 
and pab — 1 — Paa fraction of them to pre-existing nodes in B. The similar process takes 
place when a new node joins to B, where the new node has mg < m outgoing links from 
which pbb of them goes to nodes in B and the complementary pba = 1 — Pbb goes to A. 
The attachments to nodes in either networks are preferential and the rate of acquiring a 
link depends on the number of connections and the initial attractiveness of the pre-existing 
nodes. 

We define Na^a, t) as the average number of nodes with total Jva number of connections 
that includes the incoming intra-links Jvaa and the incoming cross-links ksA in network A at 
time t. Similarly, Nb^e, t) is the average number of nodes with ks = ksB + kAB connections 
at time t in network B. Notice that the indices are discriminative and the order in which 
they are used is important, as they indicate the direction that the links are made. Further 
more we also define Naa^aa, t) and NssiksB, t) the average number of nodes with kAA and 
ksB incoming intra-links to A and B respectively. Finally, we also have N BA{kBAi t) and 
Nab^ab, t) to denote the average number of nodes in A and B with ksA and /cab incoming 
cross- links. 

To keep this paper less cumbersome we will only analyse the time evolution of network 
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FIG. 1. Schematic illustration of the growth processes of network A and B. Each new node ua 
or n B is making a fractional cross linking to other network as well as intra links to their own . 

A and apply our results to network B. In addition to this, we only need to give the time 
evolution of Na^aa, k B A,t), denned as the joint distribution of intra-links and cross-links. 
Using this distribution we can find all other distributions that are mentioned earlier. The 
time evolution of Na^aa, k B A,t) can be described by a rate equation 

dN A {k A A,k B A,t) 1 { Ul 1 , , | \ at ( i ii + \ 
dt = MA~(t) Am - 1 + k BA + a)N A {k A A - 1, k B A, t) 

-(k A A + k BA + a)N A (k AA , k B A, *)] 

+PBAm B [(k A A + k BA - 1 + a)N A (k A A, k BA - 1, t) 

-(k A A + k BA + a)N A (k A A, k BA ,t)]} + h AA oh BA o- (1) 

The form of the Eq. (fj) seems very similar to the one used in [15|. In that model the rate 



of creating links depends on the out-degree of the issuing nodes and the in-degree of the 
target nodes. Here we are concerned with two different types of in-degrees namely intra- 
and cross-links of the nodes. 

On the right hand side of Eq. ([]]) the terms in first square brackets represent the increase 
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in the number of nodes with ftA links when a node with kAA — 1 intra-links acquires a 
new intra-link and if the node already has links this leads to reduction in the number. 
Similarly, for the second square brackets where the number of nodes with kA links changes 
due to the incoming cross-links. The final term accounts for the continuous addition of new 
nodes with no incoming links, each new node could be thought of as the new publication in 
a particular research discipline. The normalization factor Ma(<) sum of all degrees is defined 
as 

oo 

M A (t)= £ A kA N A (k A ,t). (2) 

k A =0 



We limit ourself to the case of preferential linear attachment rate [13 



A kA = a + k A , (3) 

shifted by a > 0, the initial attractiveness 0] of nodes in A, which ensures that there is a 
nonzero probability of any node acquiring a link. The nature of A kA lets one to obtain, as 
t — > oo 

M A (t) = (a+ < m A >)t, (4) 

where < rriA >= Paa^ia + Pba^ib is the average total in-degree in Network A. Eq. (f|) 
implying that M^(t) is linear in time. Similarly, it is easy to show that Na^aa, ksA,t) = 
^a(^4A, kBA)~t is also linear function of time. We use these relations in Eq. (JJ) to obtain the 
time independent recurrence relation 

[a+ <m A > + <m A > (k AA + k B A + a)]n A (k A A, k BA ) 
= PAAm A (k A A + k B A + a- l)n A (k A A - 1, k BA ) 
+PBAm B (k AA + k BA + a- l)n A (k A A, k BA - 1) 

+ (a+ < m A >)S kAA oS kBA o. (5) 

The expression in Eq. (^) does not simplify however, it lets us to obtain the total in-degree 
distribution 



k A 

N A (k Aj t)= N A {k AA ,k A -k AA} t). (6) 

k A A=0 

Writing N A (k A ,t) = n A (k A )t and since k A = k AA + k BA then n A (k A ) satisfies 

[a+ < m A > + < m A > (k A + a)]n A (k A ) =< m A > (k A + a - l)n A (k A - 1) 

+ (a+ < m A >)4 A o- (7) 

Solving Eq. (0) for k A > yields, 

T(k A + a)T(a + 2 



n A (0)=(l + - ^-—) . (9) 



with 



< m A > +a, 

As k A — > oo Eq. @ gives the asymptotic behaviour of the total in-degree distribution in A 

n A (k A )^k A ^, (10) 

which is a power-law form with an exponent 7a = 2 + a/ < m A > that only depends on the 
average total in-degree and the initial attractiveness of the nodes. Similarly, we can write 
the total in-degree distribution in network B for the asymptotic limit of ks as 

n B {kB) ~ ^~b b with 7s = 2 + bj < m B > . (11) 

Again, the exponent depends upon the initial attractiveness b of nodes and the average total 
incoming links < m B >= PBB m B +Vab^a- 

We now move on to analyse N AA (k AA , t), the distribution of the average number of nodes 
with k AA intra-links in network A. In citation network one can think of these links being 
issued from the same subject class as the receiving nodes and in the case of human sexual 
contact network, they represent the homosexual interactions. Since 

oo 

N AA {k AA ,t)= ]T N A (k AA ,k BA ,t), (12) 

k BA =0 



which can also be written as NAA{kAA,t) = n>AA{kAA)t, a linear function of time. Then 
summing Eq. (H) over all possible values of k B A 



nA(k A A) = n A(k A A,k B A), (13) 

k BA =0 



we get 



[a+ <m A > + <m A > (k AA + a)]n AA (kAA)+ < m A > ^ k BAn A (k A A, k BA ) 

k B A=0 

oo 

= PAAm A (k A A + a- l)n A A(k A A - 1) + Vaa^a k BA n A {kAA - 1, k BA ) 

k BA =0 

oo 

+PBAm B (k A A + a)n A A(k A A) + VBA°m B ^ (^ba - l)n A (k A A, k BA - 1) 

k BA =0 

+ (a+ < m A >)S kAA o. (14) 

For large kAA Eq. (0) reduces to 

[a+ <m A > +PAAmA(k A A + a)}n A A(k A A) = 
PAAm A (k A A + a- l)n AA (kAA ~ 1) + (a+ < m A >)h AA o- (15) 

Iterating former relation for kAA > yields 



T(k AA + a)T(a + 2 + 



p BA m B +a\ 



nAA{k A A) = ^ , , ^^0 ^(0)^ (16) 



where 



T(k AA + a + 2 + P^ m * +a )r(a) 



n^(0) = (l+ )- X . (17) 

a+ < m A > 

In the asymptotic limit as kAA —* oo Eq. (|16|) has a power-law form 

71aa(A:aa) ~ Kka aa -. witn exponent = 2 + Pba7TIb + a (is) 

that depends upon both paa and the coupling parameter p B A- 

Similarly, the time independent recurrence relation for N BB (k BB ,t) has the same form 
as Eq. flT3p with the only difference being the parameters. Therefore we will simply give the 
power-law distribution 



n BB (k BB ) ~ k BB BB with lBB = 2 + PABmA + \ (19) 

VBB-mB 

where the other coupling parameter pab is revealed in the exponent. 

Finally, the distribution of average number of nodes with incoming cross-links 
N B a (k ba 1 1~) in A can be found by summing over N A (k A A, k B A,t) for all its intra-links 

oo 

N BA (k B A,t)= £ N A (k A A,k BA ,t). (20) 

k AA =0 

As before N B A{k B A,t) = n BA (k BA )t is also linear in time. When the cross links k B A are 
large enough, then from Eq. (^) we obtain 

T(k BA + a)T(a + 2 + P*4 m *+* ) 

\ an i PBAm B ' v ' 



where 



lew ci , Pba™b& 
^ba(O) = (H ) • (22) 

In the asymptotic limit as kAB —> oo the distribution 

n BA {k B A) ~ k^ sa , 7ba = 2 H (23) 

Pbattlb 

has a power-law form and similarly for the network i? as A;^^ — > oo 

nxfl^a) ~ ^ with 7AJ? = 2 + PBBmB + 6 . (24) 

PABm A 

Unlike the case in intra-links, here the exponents are inversely proportional to the coupled 
parameters pba and pab respectively. 

III. DISCUSSION AND CONCLUSIONS 

For the sake of simplicity, we set the number of outgoing links of the new nodes in either 
networks to be the same, i.e. rriA = m B = m. Furthermore taking the rate of cross linking 
to be pab = Pba = q and the rate of intra linking paa = Pbb = P, consequently we have 
p = 1 — q, and q as the coupling parameter. 
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In the weak coupling case, the cross linking is negligibly small i.e. q — > then the 

power-law exponent of the intra-link distribution is 

qm + a . . 

lAA = 2 + f — 25 

(1 — q)m 

equal to total link distribution jaa = 2 + a/m. This gives a solution obtained in [|J and 
when a = m we recover the exponent 7 = 3, the empirical findings in 

The case of strong coupling q — > 1 is an illustrative example for citation networks, which 
results in ^aa = Ibb —> oo, both intra-link distributions having exponents approaching to 
infinity. 

Thus, varying q in (0, 1) yields any values of ^aa between (2 + a/m) and oo. On the 

contrary, the exponent of cross-link distribution 

(1 — q)m + a , . 

lBA = 2 + ^ q -> 26 

qm 

decreases from oo to 2 + a/m, as q increases from to 1. 
Taking a = m gives 

1 + 9 , 2 

Iaa = 2 + and lB A = 1 + -• 27 

1-q q 

Supposing < q < 0.5, which seems reasonable for consideration of citation networks, we 

find that 3 < ^aa < 5 and 'Jba > 5. The former result coincides with the distribution of 

connectivities for the electric power grid of Southern California Where the system 



is small and the local interactions is of importance hence there seems to be some analogy 
to the intra-linking process. For the latter, as far as we are aware there is none empirical 
studies present in the published literature. 



Now, consider the web of human sexual contacts [fTl,12l. If we let A to represent males 

and B females that is A = M and B = F then 

(1 — q)m + a , n (l — q)m + b 

lFM = 2 + ± ^ and lMF = 2 + ± ^ 28 

qm qm 

are the power-law exponents of the degree distributions of the sexes. Where a and b denote 

the male and female attractiveness respectively and usually a < b is considered [I2||. By 



setting a/m = 1.31, b/m = 1.54 and q — > 1 that is, cross links are predominant then as in 
fT2|j we obtain j FM = 1 + o>fm = 3.31 for males and ^mf = 1 + &mf = 3.54 for females. The 
exponents a™ = 2.31 and olmf = 2.54 have been observed for the cumulative distributions 
in empirical study [|TT|. 



The model we studied here seems to have the flexibility to represent variety of complex 
systems. 
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